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Low energy photon production is an important background to many current and future precision neutrino 
experiments. We present a phenomenological study of t-channel radiative corrections to neutral current neu- 
trino nucleus scattering. After introducing the relevant processes and phenomenological coupling constants, 
we will explore the derived energy and angular distributions as well as total cross-section predictions along 
their estimated uncertainties. This is supplemented throughout with comments on possible experimental signa- 
tures and implications. We conclude with a general discussion of the analysis in the context of complimentary 
methodologies. 



1. Introduction 

Recent neutrino scattering experiments report sig- 
nals with accuracies below the 1% level. Such un- 
precedented sensitivities demand corresponding ef- 
forts to determine backgrounds. Radiative correc- 
tions are clearly expected at this level. A proper 
understanding of induced photon production is espe- 
cially critical for those experiments searching for elec- 
tron neutrino appearance with non-magnetized detec 
tors Aguilar et al. 2001 , Aguilar-Arevalo et al. 2007 



2009 where it is difficult to distinguish gamma radi- 
ation from electrons. This is the case for many pre- 
cision short baseline oscillation experiments. Stan- 
dard radiative corrections from final state photon 

and resonant 



bremsstrahlung Rein and Sehgal 
A/N* production Adler 
amincd in the literature 



1968 



1981 



have already been cx- 
and are included in exper- 



imental Monte Carlo simula tions | Agostine lli et al.| 
|2003| , |Casper||2002| , [Garvey|[2009| . Next generation 
magnetized detectors will alleviate some the uncer- 



tainties induced by this background Chen et al. 2007 



In what follows, we present a novel Standard Model 
contribution to t-channel neutral current photon pro- 
duction in neutrino nucleon scattering first introduced 
by us in Jenkins and Goldman 2009 . In contrast 



to the well known s-channcl effects described above, 
our selected class of processes are less obviously con- 
nected to the external line quanta. We consider both 
neutrino and anti-neutrino processes and our results 
may be extended to other similar interactions both in 
neutral and charged current scattering. Although our 
primary focus is on modest energies, our results are 
relativistically covariant and thus may be applied to 
any energy. Of course, at high energies Regge trajec- 
tory generalizations of the meson exchanges are nec- 
essary which will naturally lead to a quark picture of 
the interaction. 

This paper is organized as follows. In section [2] we 
introduce the dominant scattering process diagrams 
and phenomenologically derived coupling constants. 
This is followed by a derivation of the scattering cross- 




Figure 1: Specific diagrams considered in this analysis. 
Variants and interference effects are discussed in the text. 



section. We show our numerical differential and total 
cross-section results in section [3] where we also point 
out the importance of interference effects. We con- 
clude in section [4] with a brief summary and general 
discussion of our methodology. 



2. Process 

2.1. Diagram and Couplings 

Figure [l] shows the two dominant t-channel modes 
considered in this analysis, differentiated by interme- 
diate u> and p° meson exchange. In both diagrams, 
the Z-boson carrying the neutral current interaction 
from the neutrino line mixes into a vector boson in the 
familiar fashion of Vector Meson Dominance Schild- 



knecht| [2006 . The hadronic vector meson then un- 
dergoes a virtual decay to a photon and pion in the 
i-channcl. This last couples strongly to the hadron 
target. Of course, there are other similar contri- 
butions from vector-meson (Regge) recurrences, but 
these predominantly affect only the overall strength 
for q 2 <g; M' 2 , where the excited state is integrated 
out of the interaction. Low energy hadron scattering 
experiments suggest that at modest energies the sum 



2 
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over all such contributions is likely to be dominated 
by these leading ones. For the remainder of this sec- 
tion we will focus on the the to exchange diagram and 
discuss the effects of interference in subsection 13.21 

The anatomy of this diagram is shown in figure [2j 
where the needed coupling constants are circled for 
emphasis. These are extracted phenomenologically 
from measured processes. 

Beginning with the 7r — 7— meson vertex. We 
see that this contribution is similar to the triangle 



anomaly mediated interaction identified in Harvey 
et al.| [2007] and discussed in [Hiil] |2009] . Our advan- 
tage over this approach is that the vertex strengths are 
known phenomenologically from the decay lu — > it + 7 
computed from the effective Lagrangian term 



eg, 



(1) 



Here is the photon field strength tensor and the 
electromagnetic coupling e is factored out for conve- 
nience since it is necessarily present from the photon 
interaction. We point out that, although Eq. has 
the same form as that induced by the triangle anomaly 
due to the axial vector nature of the pion current, it 
exists independent of the anomaly. 

Using this interaction, and neglecting the ir° mass, 
we find the squared decay amplitude 



A 2 = - 



~^~ k ' q ~ 6 ■ 



(2) 



where k and q are the photon and pion momenta, 
respectively. This implies the decay width 



T(UJ — > 7T + 7) 



24 



(3) 



Fitting Eq . 
eTaD [27)06], 

L87M 



(|3j to the observed decay width Yao 
we extract the coupling constant g W77r — 





Figure 2: Phenomenologically extracted coupling con- 
stants determined from experimental data. 



find 



9u,z{v 2 ) 



\ / dz 

'0 -In 



= -gg^ q gM 2 s 2 v 

I2ir 2 cw 

1 ' p 2 z(z — 1) + TIT, 2 



(4) 



dx 



p 2 z(z — 1) + m 2 + x(M 2 — m 2 ) 



^similar exercise may be performed with the 
p decay, in which case one extracts g PT n — 0.55 /M p . 

The measured partial decay widths allow for very 
accurate coupling constant extraction beyond the 
0(10%) level shown here. For the purposes of describ- 
ing a sub-1% signal, our accuracy adequately provides 
total cross-section predictions to better than 0.1%. 
This reasoning holds for other parameter extractions 
given throughout the text. 

Moving on, the strength of the Z — lu mixing and 
its p 2 running may be extracted from the self energy 2 2 CfOSS SGCtiOfl 



after dropping logarithmically divergent contribu- 
tions. Here cw and sw are the cosine and sine of 
the weak mixing angle. Taking reasonable limits of 
Eq. p| yi elds simplified analytic results Jenkins and 
I Goldman] |2009] but the remaining Feynman integrals 
may be easily performed numerically. Doing this, 
we find the averaged g^z = 600 MeV 2 assuming 
m q ~ 3MeV and M ~ M u at momenta transfers be- 
tween 200 - 1000 MeV. We find a slight O(10%) vari- 
ation of g^zip 2 ) within this region of interest. 

For the remaining couplings, we make use of the 
Standard Model weak interaction of the Z-boson to 
neutrinos and quarks and the well known pion cou- 



pling to the nucleon Liu et al.| 1995 via the interaction 
Ci = g 7r jv7v^ , 7 M 7 5 i9 At 7r-T^', where \t is the nucleon field 
and t are isospin generators. No other parameters are 
required, so the prediction of the contribution to the 
total cross-section for producing a final state photon 
is absolute for this diagram. The analogous analysis 
is easy to perform for the p° exchange case. 



diagram shown in figure [2] Following Goldman et al. 
1992ajb| we parameterize the uj — q — q form factor by 
9uqqM / (M 2 —p 2 ) where M describes the finite size of 
the uj meson. The bare uj—q—q coupling is fo und to be 
g^qq ~ 3.1 from the decay ui 



Evaluating the u> exchange diagram of figure [TJ we 
find the squared scattering amplitude 



Yao et al. 



2006 . Calculating the self energy via dimensional 
regularization and considering only those terms that 
contribute to the p 2 dependence of lu — Z mixing we 



A 2 



l28M 2 N gl z g 2 gl NN gl z (q 2 ) 



(q 2 ~ M 2 ) 2 (q 2 ~ M 2 ) 2 (Q 2 - M 2 ) 2 
x ti ■ t f {pi ■ Pf - M 2 N ) ((k ■ i,) 2 + (k ■ £ f ) 2 ) 



(5) 
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in terms of the labeled four momenta. The upper 
portion of Eq. ^ shows the general coupling con- 
stant and propagator dependencies while the lower 
factor describes the kinematics that follow from the 
diagram's Lorentz structure. In the Center of Mass 
(CM) frame, the momenta can be written explicitly 
as 



Pi 

If 
Pf 
k 



{E Pi ,—Ei i ) 
{E lf , Ei f ) 
{E P f , p Pf ) 

(E k , Ek), 



(6) 

(7) 

(8) 
(9) 

(10) 



where we employ the shorthand Ei — pi to indicate 
a massless particle's 3-momentum of magnitude -E 1 ,. 
In this frame we find, after performing trivial integra- 
tions over momentum-conserving delta functions, the 
differential cross-section for the final state photon's 
energy and angular distribution to be 



do 



dEkd[i 
dEf f < 



{2tt) 4 E u (E u +E Pi ) 
9lz(q 2 )q 2 Q 2 ((* • h 



(11) 



(k • if) 



( q 2 - M 2 ) 2 ( q 2 - m 2 ) 2 (q 2 - M„y 



where the momenta transfers are given by q 2 = 
-2E u E lf {l - fj,i f ) and Q 2 = q 2 - 2k ■ i t + 2k ■ £ f . 
Here [i — cos is the photon opening angle from the 
beam direction and \xi f is the cosine of the opening an- 
gle between the neutrino in the final and initial state. 
It is related to \i and the cosine of the opening angle 
between the photon and the final state neutrino (Xi f k 

by 



tilt 



1 



■ f$ k cos (j>. (12) 



This is the only <f> dependent term in the system. Mo- 
mentum conservation then fixes the remaining open- 
ing angle to be 



(i>e f k 



1 

2Ei } Ek 



(13) 



(a - 2y/i(Et f + E k ) 



2EkEe f — M N ) 



where s = {Ei i + E p .) 2 is the relativistically invariant 
squared CM energy. Additional constraints and limits 
of integration are found by requiring that fJ,i f k and Ek 
take on physical values. 



Evaluating Eq. (Ill subject to these constraints in 
the reasonable limit \q 2 \ < M§, \Q 2 \ > M 2 and 
9ujz(q 2 ) ~ 9ljZ, we integrate over (f> to obtain 



do- _ M N E2 k9lz9tz9t 1 ^9iNN 
dE k dfJ, ~ (2ir) 3 E h (E u + E V% )M% 



(14) 




Figure 3: Angular and energy differential cross-section dis- 
tributions in the CM and lab frames for various neutrino 
beam energies. 



dE l} (e 2 . (1-h) 2 + E 2 f (l - M , /fc ) 2 ) 

1 f a — b c 3 — 2ac 2 + abc — b + a 



f 2 (b~c) 2 I (62 _ i)i 
where 



(C 2 -1)« 



/ = 



2Et i Et f y/l-n*Jl-ni fl 
2E e .Ei f (l - (ifitfk) 



(15) 
(16) 



2E k Ei i (l - M ) - 2E k Ei f (l - w fc ) 
a H — (17) 



Ml 
f 



(18) 



Assuming physical parameters, the dimensionless 
quantities a, b and c arc all greater than unity. This 
leaves only the one-dimensional integral over the final 
neutrino energy (EiA to perform. 



3. Phenomenology 



In what follows, we numerically explore Eq. (Ill in 



both the CM and lab frames using the phenomcnolog- 
ically derived coupling constants. We point out that 
we are using the full cross-section expression without 
approximation, including the q 2 running of g u z{q 2 )- 
We first discuss the results of the oj exchange process 
alone followed by the modifications induced by inter- 
ference with the contribution of the p° graph.. 



3.1. Results 

In the CM frame, the predicted cross-section is 
weakly peaked in the backward direction with an en- 
ergy maximum near the highest kinematically allowed 
energies due to the overall factor of E k in Eq. ( |14[). 
This can be seen in the upper panels of figure I3J 
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Figure 4: Lab frame differential cross-section contour plot 
with beam energy Ee i = 200 MeV. Angular and energetic 
projections are shown for convenience. 



350 MeV 




Figure 5: Lab frame differential cross-section contour plot 
with beam energy Ei i = 350 MeV. Angular and energetic 
projections are shown for convenience. 



Figure 6: Lab frame differential cross-section contour plot 
with beam energy Ei i — 500 MeV. Angular and energetic 
projections are shown for convenience. 



1000 MeV 




Figure 7: Lab frame differential cross-section contour plot 
with beam energy Ee i — 1 000 MeV. Angular and energetic 
projections are shown for convenience. 



Boosting these distributions to the lab frame pushes 
the angular distribution forward and spreads the en- 
ergy of the photon, as is evident in the lower pan- 
els. Numerically integrating Eq. (14 1, we plot the 



still affect current Aguilar-Arevalo et al. 2007 2009 



lab frame differential cross-section for beam energies 
of 200 MeV, 350 MeV, 500 MeV and 1000 MeV in 
figures [4] [5] [6] and [7] respectively. In each case, we 
display the E 1 and cos 9 dependent contour plots as 
well as energy and angular projection panels obtained 
by integrating over one of the variables. The total 
cross-section is also noted for reference. The angu- 
lar distribution moves toward the forward peak with 
increasing neutrino energy due to the growing boosts 
from the CM to the lab frame. The distribution con- 
sistently peaks near the center of the kinematically 
allowed photon energy range. 

Integrating over the final state photon energy and 
angular distribution, we plot the total cross-section 
as a function of neutrino beam energy in 

figure m At 

high energies the cross-section grows as \JE V and near 
threshold as E%. A logarithmic insert plot showing 
the low energy region of interest is included for conve- 
nience. Here, the cross-section is roughly three orders 
of magnitude smaller than the typical charged current 
cross-sections Lipari et al. 1995 . However, this may 



and future Chen et al. 2007 experiments. Addition- 



ally, long baseline and precision scattering neutrino 
experiments performed at higher energies (see, for ex- 



ample 


Adams et al. 2008 , 


Adamsonet al. [2008|, 


Aoki 


et al. 


2003 , 


Ayres et al. 


2004 , 


Beavis et al. 


2002 



and references therein) will be sensitive to this class of 
processes with an order of magnitude enhanced cross- 
section. 



3.2. Interference Effects 

We now discuss the interference effects resulting 
from the addition of the p° exchange diagram of figure 
[l] We point out that the p° exchange mode's scatter- 
ing amplitude will have the same form as in the u) case 
with different (but still phcnomcnologically extracted) 
coupling constants and exchanged masses <-> M p o . 
The couplings only effect the overall magnitude while 
the meson masses can also influence the cross-section 
shape. Since M L 



M„ 



2006 



we see that 
d be almost 



Yao et al 

the resulting spectral distributions shou 
identical. Thus, it is sufficient to consider variations 
between the overall scattering magnitudes induced by 
coupling constant differences. 
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Figure 8: Total cross-section as a function of neutrino 
beam energy. A log scale insert plot is given to emphasize 
the low energy region of interest. 



Such relative differences occur due to the 7r — 
7— meson vertex as well as in the meson-Z mixing 
term. In the first case, the relevant coupling constants 
were calculated in subsection |2. 1 1 leading to a suppres- 
sion 



9p~fx 



0.55M, 



g UJ7r 1.8M p 



= 0.31. 



(19) 



The meson- Z mixing contribution is less trivial to 
understand. On the basis of SU (3) flavor symmetry, 
one expects similar results for the p° — Z and u> — Z 
mixing terms up to isospin effects. Looking at the self 
energy diagram in figure [2] we see that the u> couples 
equally to the u and d quarks that contribute to the 
loop, whereas the p° does so with opposite signs due 
to isospin. The standard model couplings are 



gz-u 



9zdd 



9 

4:C W v 3 
Ac w 3 w ' 



(20) 
(21) 



As might be expected from the fact that the Z-boson 
is dominantly isospin one like the p, the Z — p mixing 
is enhanced relative to the Z — lu mixing. Combining 
this reasoning with SU(3) flavor symmetry breaking 
manifest in deviations of g pq ql g^qq from unity, we find 




gpqq 

g^jqq 



2s 2 



tttt) <j>{u) — > 7T7r7r) 
7T7r7r) <j)(p — > 7T7r) 



(22) 



= 4.1, 



where the </>'s denote phase space integrals required 
for the meson— q — q coupling constant extractions. 
Thus, we estimate comparable cross-sections given 

by 

a p /a u ps (5p- 7 - 7r /.g w -7-7r) 2 x (g pZ /guz) 2 ~ 1-6, 

(23) 



which arises from the accidental cancellation of the 
meson— 7r — 7 suppression and the meson-Z mixing 
enhancement. 

The amplitudes for these processes are similar and 
significant interference is expected to occur. From 
this effect the overall cross-sections may be modified 
by a factor between 0.07 and 5.1 for total destructive 
and constructive interference respectively. Within the 



framework of the triangle anomaly Harvey et al. 



2007 



one may gain a handle on the relative interference 
phase by considering the low energy limit where the 
u> and p contributions must sum to yield the ir° — 7 — Z 
coupling equal to 1 — 4s^, . This is small and picks out 
the lower bound of our interference region. However, 
since the phase relations of our phcnomenological am- 
plitudes are not fixed and are independent of the tri- 
angle anomaly, isospin constraints from the quark cou- 
plings to the Z-boson cannot be applied. This allows 
for additional latitude in matching experimental re- 
sults. Cross-sections yielded by the lower limit fall 
well below expected future experimental sensitivities 
and therefore predict a negligible background. Con- 
tributions at the upper limit would have a substantial 
observable impact on precision measurements, and as 
such, should be included in future experimental Monte 
Carlos. 



4. Conclusions and Outlook 

Other processes related to those in figure [l] such as 
by the exchange of the to and 7r° for other mesons 
(with the same quantum numbers), will contribute 
to similar production of photons in the final state. 
These will differ from our calculation only by the cou- 
pling constants and meson masses which are neces- 
sarily heavier and should lead to propagator suppres- 
sions. 

After studying the diagrams in figure [l] it is clear 
that this process class yields identical results for both 
neutrinos and anti-neutrinos. The only difference be- 
tween these amplitudes resides in a sign change at 
the axial-vector neutrino-Z coupling. This vanishes 
when Lorentz contracted with the rest of the diagram 
which is symmetric under the free indices at the ver- 
tex. Additionally, we find by means of direct compu- 
tation that many variants of figure [l] vanish due to 
similar symmetry reasons. In particular, amplitudes 
from diagrams with axial vector or pseudoscalar, as 
opposed to vector, meson exchange vanish. Addition- 
ally, the "reversed" diagram, where the ir° couples to 
the neutrino line, yields a null contribution. We point 
out that such a contribution, if nonzero in principle, 
would be highly suppressed by the tt — v — v coupling. 

Throughout this analysis we have used an on-shell 
coupling strength for the u) — tt — 7 vertex, which is 
a commonly used phcnomenological technique. The 



G 
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slow variation of the computed Z — u) mixing sup- 
ports such an approach, but a three body interaction 
could behave differently. In this case one would ex- 
pect a decrease in amplitude as vertex form factors 



act to suppress the effective coupling Lansberg 2007 
The vertex structure for the to — q — q coupling, also 
obtained from an on-shell decay, produces the same 
kind of uncertainties. We believe that this issue is not 
a serious problem for our analysis, as the Z — lu mixing 
is only used within a few mass squared units from the 
on-shell point. 

Another potential background in is e appearance 
searches occurs when the a decay photon from a pro- 
duced 7T° is lost to the detector, leaving a single photon 
faking an electron track. Fortunately, in this case the 
event rate can be normalized to the corresponding pro- 
cess in charged current neutrino scattering, which pro- 
duces a neutral pion in conjunction with the charged 
lepton. Although this may dominantly occur due to 
intermediate state processes, such as production of a 
A baryon or N* followed by its decay back to a nu- 
cleoli and a pion (see Paschos et al. 2000] and ref- 
erences therein), concern also arises regarding other 
processes, including those that may be coherent over 



the entire nuclear target with an amplified rate Isik- 
sal et al.| |1984] , |Rein and Sehgal| |1983] . The pre- 



ceding analysis may be applied, in an analogous way, 
to coherent f-channel pion production. This parallel 
process is interesting from the point of view of inter- 
ference. If the interference is destructive for photon 
production, there can be a significant difference be- 
tween coherent pion production between charged and 
neutral current modes, as this interference cannot oc- 
cur in the charged current case. We will explore this 
possibility in a future study. 
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